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ABSTRACT 

We explore the dependence of the central logarithmic slope of dark matter halo density 
profiles a on the spectral index n of the linear matter power spectrum P(k) using cosmological 
TV-body simulations of scale-free models (i.e. P(k) cx k n ). These simulations are based on 
a set of clear, reproducible and physically motivated criteria that fix the appropriate starting 
and stopping times for runs, and allow one to compare haloes across models with different 
spectral indices and mass resolutions. For each of our simulations we identify samples of well 
resolved haloes in dynamical equilibrium and we analyse their mass profiles. By parameterising 
the mass profile using a "generalised" Navarro, Frenk & White profile in which the central 
logarithmic slope a is allowed to vary while preserving the r~ 3 asymptotic form at large 
radii, we obtain preferred central slopes for haloes in each of our models. There is a strong 
correlation between a and n, such that a becomes shallower as n becomes steeper. However, if 
we normalise our mass profiles by r_2, the radius at which the logarithmic slope of the density 
profile is —2, we find that these differences are no longer present. This is apparent if we plot 
the maximum slope 7 = 3(1 — p/p) as a function of r/r_2 - we find that the profiles are 
similar for haloes forming in different n models. This reflects the importance of concentration, 
and reveals that the concentrations of haloes forming in steep-n cosmologies tend to be smaller 
than those of haloes forming in shallow-n cosmologies. We conclude that there is no evidence 
for convergence to a unique central asymptotic slope, at least on the scales that we can resolve. 
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1 INTRODUCTION 

The currently favoured model of cosmological structure formation 
asserts that the Universe is dominated by some form of non-baryonic 
Cold Dark Matter (hereafter CDM), and it makes a number of 
generic predictions about how the dark matter clusters on small 
scales. Arguably the defining prediction of the CDM model is that 
the radial mass density profile of dark matter haloes is divergent at 
small radii, which leads to a central density "cusp". Characterising 
the form of this mass density profile has been one of the most active 
research problems in computational cosmology over the last decade. 

The seminal study of Navarro, Frenk & White (1996,1997; 
hereafter NFW) introduced the concept of a "universal" mass density 
profile for dark matter haloes that form in hierarchical clustering 
cosmologies. This NFW profile is written as 



p(r) = 



p(r) 



(r/r 3 )(l + r/r 3 y 



(1) 



where r s is a scale radius and p c is a characteristic density, which 
can be related to r s once the virial mass of the halo is fixed; equa- 
tion[TJdescribes a one-parameter family of curves. It is convenient 
to rewrite equation[T]in the more general form 



(r/r s ) a (l + r/r s ) 3 -" 



(2) 



where a is the central asymptotic logarithmic slope. The NFW pro- 
file is "universal" in the sense that it describes the ensemble averaged 
mass profile of dark matter haloes in dynamical equilibrium, inde- 
pendent of virial mass, cosmological parameters and initial power 
spectrum. 

During the last decade numerous studies have investigated 
whether or not the NFW profile does indeed provide an adequate 
description of the mass profile of dark matter haloes, and to under- 
stand the physical mechanisms that shape the functional form of the 
profile. NFW recognised that r s correlates with the virial mass of 
the halo, increasing with decreasing virial mass. They argued that 
this virial mass-scale radius relation is an imprint of the hierarchical 
assembly of haloes; low-mass haloes tend to collapse before high- 
mass haloes, when the mean density of the Universe is higher, and 
r 3 reflects the mean density of the Universe at the time of collapse. 
The halo sample that formed the basis of |Navarro et al.| 



( 1997 1 contained of order 10 particles (and hence resolving the 



profile down to about 10% of the virial radius, according to the 
convergence criteria of P ower et al.| (j2003 >). Subsequent studies 
drew upon haloes containing about two orders of magnitude more 
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particles within the virial radius and confirmed the basic finding 
of NFW that CDM haloes are cuspy (e.g. Fukushige & Makino 
1997l|Moore et al.|1998||jing & Suto|2000||Jing|2000||Fukushige| 
& Makino|2001||Jing & Suto|2002fr. These studies led to debate 



about the exact value of the logarithmic inner slope, which ranged 
from a ~ 1 to 1.5 (e.g., |Navarro et al.|199"7[|Moore et al.|1999[ >. 
The study of Power et al. ( 2003 1 provided convergence criteria that 
allowed simulators to understand the impact of numerical artifacts 
on the central structure of haloes, and to identify the innermost 
radius at which the mass profile could be considered reliably 
resolved. More recently, a great deal of effort has gone into even 
higher resolution simulations that have revealed density profiles 
of CDM haloes to well within 1% of the virial radius (Fukushige 
et al.|2004HTasitsiomi et al.|2004| [Navarro et al.|20041|Reed et al.| 



2005 |Diemand~et al. 2005 ). The highest resolved simulation to-date 
reached an effective mass resolution of about 130 million particles 
in a cluster sized dark matter halo still supporting the evidence for 
a central cusp with a logarithmic inner slope of about a ~ 1.2 
Piemand et al.|2005) . 

Cosmological simulations allow us to characterise the func- 
tional form of the density profile and to explore the important physi- 
cal processes (e.g. merging, smooth accretion) that drive this form 
(cf., |Lu et al.|2006l|Salvador-Sole et al.|2007) , but a theoretical un- 
derstanding of the origin of the mass profile is essential. Is the form 
of the profile set by non-linear processes during the virialisation of 
the halo, or is there an imprint of the primordial power spectrum 
P(k)7 

One must make strong assumptions about the connection be- 
tween density and velocity dispersion to obtain analytical predictions 
about halo structure. For example, under the assumption that the 
phase-space density is a power law in radius as suggested by |Taylor| 
|& NavarrolpOOT) one can solve the spherical Jeans equation to 
obtain the density profile ([Hansen 2004; Dehnen & McLaughlin 
|2005| [Austin et al.|2005[ >. These studies, guided by the results of 
the numerical simulations, confirm the central logarithmic slope to 
be in the range a ~ 1 to 2, although Hansen & Stadel (2006) have 
claimed that equilibrated haloes have central slopes of a ~ 0.8. 

However, it is interesting to ask whether or not there is a de- 
pendence of halo structure on the primordial power spectrum P(k). 
In particular, does the form of the profile depend on the (effective) 
spectral index n of P(fc)? One of the key predictions of the NFW 
papers is that the shape of the "universal" profile should hold in 
any hierarchical cosmology, including scale-free ones. This was 
confirmed by subsequent numerical simulations (e.g. Cole & Lacey| 
|1996[ >. However, a number of analytical studies have claimed that 
the density profile should depend on spectral index (e.g., Hoffman 
& Shaham|T985] |Syer & White|1998l|Subramanian et al.||2000f 
Huffenberger & Selja k|2003[|Salvador-Sole et aLpOOTp ind these 



have been supported by numerical simulations providing similar 
evidence (e .g., |Crone et al.|1994[|Eke et al.|200l||Ricotti|2003"][Cenl 



|et al.|2004||Reed et al.|2005[|Ricotti et al.|2007| >. 

Within this context it is interesting to consider the findings 
of |Ricotti| ( [20031 and |Cen et "at] pOOl} . Both studies analysed 
high-resolution cosmological simulations and each argued that dark 
matter haloes at high redshifts have shallow central logarithmic 
slopes, with a ~ 0.2 — 0.4. In particular, Ricotti (2003 ) claimed 
that the central logarithmic slope depends explicitly on the effective 
spectral index n, in agreement with the predictions of |Subramanian| 
|et al.| ( |2000| >. We note that |Colm et al] ( |2004| > performed and 
analysed simulations comparable to those presented in |Ricotti| 
( |2003| > and found no evidence for the shallow "cores" (a ~ 0.2); 



rather they found that their data favoured "cusps" (a w 1). 
However, Ricot ti et al.| ( |200"7] > have recently revisited this topic 
using simulations with higher mass and force resolution, and they 
continue to argue in favour of shallow cores, in agreement with 
their earlier study. 



[Ricotti (2003]) and |Ricotti et al.| l |200"7] > sought to isolate the 
effect of spectral index n on halo structure by using a series of sim- 
ulations of the ACDM model with different box-sizes and evolved 
for different numbers of expansion factors, to mimic a single n. In 
this short paper we revisit this topic and the claims of Ricotti ( 2003 1 
and Ri cotti et al.| |2007 ) using high resolution cosmological TV-body 
simulations of scale-free models in which we vary systematically the 
spectral index n. This "clean" approach allows us to study whether 
the density profiles of dark matter haloes indeed show a dependence 
on the spectral index n. 

In §[2] we give a detailed description of how we have set up and 
run our scale-free simulations. There are some important issues to be 
considered when deciding on when to start a scale-free simulation 
and when to stop the simulation to compare halo properties. We 
encapsulate our findings in a set of well defined and physically 
motivated criteria. In §[5] we present the results of our analysis of 
the halo mass profiles, and we summarise our results in § |4] 



2 NUMERICAL SIMULATIONS 

As we intend to perform scale-free cosmological simulations, we 
use a power law for the initial power spectrum of the simulation 

k 



P(k) = a AP* 



Ny 



(3) 



with A being the amplitude scaling factor and 
Avn = ( t= J and fc Ny 



Ntv 



denoting the power of the white noise of a random distribution 
of N 3 particles in a box with size L and the Nyquist frequency, 
respectively. The simple scaling with expansion factoQ a is valid 
only for an Einstein-de Sitter universe, which is the cosmology of 
our choice. The box-size L is completely arbitrary because of the 
scale-free nature of the power spectra, and so we set it to 1 Mpc. 
N denotes the number of particles along one dimension. We run 
a sequence of si mulations with the parallel A-body code Gadget2 
( Springel||2005 1 containing 512 3 particles that differ only in the 
value of the spectral index n; we use n=-0.50, -1.50, -2.25, -2.50 
and -2.75. 

In order to set-up and run the simulations, we still need to 
specify the amplitude scaling factor A for all models and to find a 
criterion for when to stop and compare objects from different runs. 
These two choices will be presented and discussed in the following 
two subsections. 



2.1 Initial Conditions 

For setting up the initial conditions of the simulations, only the 
amplitude scaling factor A in equation [3] needs to be determined; 
all other factors are given by the specifics of the model under in- 
vestigation, i.e. N and n. The mass fluctuations Ob ox within the 
computational domain can be written as follows 



We chose a = 1 at the beginning of the simulations. 
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Table 1. The parameters for the initial conditions for the simulations. 



Figure 1. Shown are the theoretical power spectra and their respective 
realisations in the initial conditions of the simulations we performed. 



2 

<?"bo 



dfe P (k) 



(4) 



where the integration range in k contains all significant frequencies 
reproduced in the box, e.g.. k x . y . z £ [2n/L, Nn/L]. With the 
definition of P(k) (cf. eq.[5} this can be evaluated analytically to be 



A = 



27<- 2 Obox n + 3 

^Ny-PWN 



(*)' 



(5) 



Here we use V2 rather than V3 in the denominator because surface 
diagonals occur more frequently than volume diagonals. Hence by 
specifying Cbox we also fix the amplitude scaling factor A. For our 
simulations we choose <7box to be 0.15 and therefore start all models 
with the same "integral power", independent of spectral index n. 

This choice for the normalisation of the initial power spectrum 
deviates from the usual course of action to fix the power at the 
Nyquist frequency &N y at the white noise level. Our experience tells 
us that such runs start at too late a stage, especially for P(k) with 
a rather steep n, because the mass variance (Jbox can be large and 
possibly non-linear. We note that previous studies have arbitrarily 
lowered the amplitude to 1/4 of the white noise level for such 
models to circumvent this problem (Efstathiou et al. 1988 , Lacey & 
|Cole|1994||Navarro et al.|1997> . 

Table[T]lists the amplitudes derived from equation[5]in terms 
of the white noise level Pwn- Figure[T|verifies that we recover the 
input power spectra from our initial conditions for the selection of 
the models under investigation. 

2.2 Stopping Criterion 

In order to define an end-point for our scale-free simulation that 
allows for a cross-comparison of identified haloes amongst runs 
with different n values we follow Na varro et al.|(T997t and track 
the evolution of M*(a), the typical collapsing mass of an object in 
our simulation at expansion factor a. This "non-linear" mass scale 
is defined by requiring that the variance of the linear overdensity 



model 


n 


A 


«r 


512-0.50 


-0.50 


0.0358 


826.0 


512-1.50 


-1.50 


0.0215 


161.5 


512-2.25 


-2.25 


0.0109 


43.9 


512-2.50 


-2.50 


0.0076 


27.3 


512-2.75 


-2.75 


0.0046 


16.7 



A is the amplitude and n the slope of the initial power spectrum. The quantity 
Of gives the adopted expansion of the universes that will lead to about 42000 
particles in an M* halo. 

field, smoothed with a top-hat filter enclosing a mass M = M*(a), 
should equal the square of the critical density threshold for spherical 
collapse (e.g., |Press & Schechter|1974||Navarro et al.|1997} 

"V. = Serif (6) 

The mass variance is related to the matter power spectrum 
P(k) via 



2 



I 

2^ 



P(k)W 2 (kR)k 2 dk, 



(7) 



where W(x) = ^(sinx — xcosx) is the Fourier-transform of the 
top-hat window function, and R — (3M/47rp cr i t ) 1//3 relates the 
spatial scale to a mass scale. 

M, (a) is a monotonic function of the cosmic expansion factor 
a and has a different shape depending on the choice of the spectral 
index n. We still need to select a fiducial value M»,f (independent 
of n) defining at[n) via inverting M*^ = M* (at{n)) for each 
of our scale-free models. A careful investigation of the maximal 
allowed expansion factor in the n = —2.5 model leaves us with 
M*,f « 42000 particles (please refer to Appendi x |Al for more 
details). The adopted at (n) values are listed in Table \\y 



3 ANALYSIS 

In this section we first discuss the tool employed to identify haloes 
(section [3~T] l and the selection criteria used to define the sample of 
virialised haloes (section [3^2] ). We then employ two approaches to 
characterise halo structure - a straightforward parameter fit to the 
density profile to assess the dependence of a on n (section |3~3) , and 
a non-parametric estimate of central concentration of the halo based 
on r_2, the radius at which the logarithmic slope of the density 
profile is —2 (section [34) . These measures of the mass profile allow 
us to explore the dependence of inner slope on the spectral index. 



3.1 Halo identification 

We used an MPI version of the AMIGA Halo Findei[^](AHF, successor 
of MHF introduced in |Gill et al.H2004"| >) for identifying haloes and 
computing their integral and radial properties. AHF locates haloes 



2 Please note that a,f for the n = —2.75 model is slightly larger than the 
maximal allowed expansion for that model as we used the n = —2.5 model 
as the reference for M„ f (cf. Appendix |A|. We therefore do not consider 
the results derived from this simulation trust-worthy but chose to present 
them anyways. 

3 AMIGA is freely available from download at |http : / /www . aip ■ de7| 
[People/AKnebe /AMIGA/] 



4 S.R. Knollmann, C. Power and A. Knebe 



Table 2. Definition of various halo samples. 



Run 


all 


M- t sample 


high-mass sample 


512-0.50 


276 


186 


38 


512-1.50 


172 


107 


44 


512-2.25 


119 


62 


33 


512-2.50 


64 


35 


19 


512-2.75 


38 


15 


16 



The columns give the number of haloes that satisfy our selection criteria. 
The first column gives the total number of objects consisting of more than 
3.15 X 10 4 particles; we refer to it as "all". The second column gives the 
number of objects in our M„ sample, i.e. haloes with masses in the range 
M £ [0.75, 1.5]M*. Finally we list the number of haloes with more than 
N > 10 5 particles, which we refer to as "high-mass sample". Note that 
all haloes satisfy our the centre-of-mass offset criterion and are therefore 
considered to be in dynamical equilibrium. 

as peaks in an adaptively smoothed density field using a hierarchy 
of grids and a refinement criterion that is comparable to the force 
resolution of the simulation (i.e. 5 particles per cell). Local potential 
minima are calculated for each of these peaks and the set of particles 
that are gravitationally bound to the peaks are identified as the 
groups that form our halo catalogue. Each halo in the catalogue 
is then processed, producing a range of structural and kinematic 
information. Haloes are defined such that the virial mass is M v i r = 
47Tp cr it Avirf5i r /3, where the overdensity criterion A v i r = 178 for 
an Einstein-de Sitter cosmology. 



3.2 Halo selection 

We restrict our analysis to only those haloes that contain in excess 
of 3.15 x 10 4 particles(~ 0.75A4,) within the virial radius. Fur- 
thermore we employ a criterion to exclude objects that we do not 
expect to be in dynamical equilibrium. For this we use A r , the 
displacement of the centre of mass of all material within the virial 
radius with respect to the centre of potential of the halo, normalised 
to the virial radius of the halo; 



A r = 



(8) 



We reject all haloes for which A r > 0.05. This is a more conser- 
vative requirement than has been used in other recent studies (e.g. 
|Neto et aTp 007 1. However, we show in another paper (Power et al, 
in prep.) that this ensures that a halo is unlikely to have experienced 
a merger with mass ratio greater than 10% in the previous dynamical 
time (defined at the virial radius). 

Additionally we define a subset of the catalog, consisting of 
haloes at about M*: we use objects in the mass range of 0.75A/* 
to 1.5M, which corresponds to halos with a number of particles 
N G [3.15, 6.30] x 10 4 . In table[2]we give a summary of the number 
of haloes used in our analysis. 



3.3 The Preferred Slope a 

We begin our study of the dependence of the central logarithmic 
slope a on the spectral index n by fitting the generalised NFW 
profile (cf. equation[2]( to each halo in our A/« sample. Recall that 

Pc 



p(r) = 



U+i 1 



(9) 
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Figure 2. The median of the quality in the different simulation boxes for 
objects in the Af* sample. 



where r 3 and p c are the fitted parameters. We perform several fits 
for each profile, varying a in the range of a — 0.3, . . . , 1.8. The 
case of a — 1.0 corresponds to the NFW profile. 

For each halo we determine a set of parameters for each choice 
of a. By comparing the \ 2 -value we can, in principle, find the 
preferred a for a halo. However, we are interested in whether or not 
a depends on n for a typical halo, and so a statistical analysis is 
required. Therefore, we define a quality measure of the fit, 



El ^ij-pfit(rj) 
bins ^ Xjj 

dof 



(10) 



here dof denotes the number of degrees of freedom of the fit, Xy 
the density value of the jth halo in the ith radial bin and p& t the 
corresponding value of the analytic function at radius n . 

Taking the median Q in each simulation box for each value of 
a we can explore how the preferred inner slope of the density profile 
depends on the spectral index n, as shown in Figure|2] Each curve 
indicates how the quality of the fit varies with a for a given n, and 
the minima of the curves pick out the preferred a for a particular 
n. For the n — —0.5 run we find a minimum at a ~ 1.6, while 
for the n = —2.75, the preferred a ~ 1.2. Even though the curves 
for the n = —1.5 (—2.25, —2.5) runs seem to have their minima 
at a ~ 1.4, a trend can be seen. Therefore we might conclude that 
the preferred a decreases with decreasing n; that is, the steeper the 
spectral index, the shallower the central logarithmic slope. 

However, care must be taken when interpreting slopes. Close 
inspection of equation|9]will reveal that the logarithmic slope will 
depend on the scale radius r s . Therefore, an apparent trend in a with 
n may instead be attributable to a trend in r s with n. If r s is system- 
atically smaller in a model with spectral index m when compared to 
a model with spectral index rii, haloes will be systematically more 
concentrated and so fits to the resolved part of the density profile 
will produce higher "effective" slopes a. We investigate this further 
inS,[34l 
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3.4 The Maximum Slope 7 

A useful non-parametric measure of halo structure is the maximum 
asymptotic slope, 



T (r) = 3 1 



p(r) 



(11) 



which uses the local mass density p and the enclosed mass density 
p to derive an upper limit to the logarithmic slope of the density 
profile at radius r. Equation| 1 1 |assumes that the halo is spherically 
symmetric with a density profile interior to r given by a power-law, 
p(r) oc r -7 . 7 defines an upper limit to the slope; a steeper slope 
would require more mass interior to r than is measured. We note 
that this measure was used by N avarro et al.| ( |2004) for resimulated 
haloes of different masses but comparable particle resolution in their 
study of the universality of the mass profile. 

In Figure [3] we plot the radial variation of 7 for all haloes 
that satisfy the selection criteria of § |3.2| and are in the same mass- 
range as those used in section [3~3| colour-coded according to the 
spectral index n of the model in which they form. Note that we have 
normalised these profiles to r_2, the radius at which the differential 
mass profile pr 2 reaches a maximum. For a NFW profile, r_2 
is identical to the scale radius r s and so it provides an attractive 
non-parametric measure of concentratioij^] 

When normalising the radius in this manner, we find excellent 
agreement between the average shapes of the 7 profiles between the 
different n models. However, the scatter between profiles within a 
given simulation is significant (cf. upper panel of figure[3]l, which 
strengthens our argument that it is essential to use a statistical 
sample of haloes when discussing the asymptotic inner slope. It is 
noticeable that the average profile in each model we have looked at 
continues to becomes shallower with decreasing radius, without 
showing evidence for convergence to an asymptotic value (cf. 
|Navarro et al.|2004] >. We find similar behaviour when considering 
only haloes in the high-mass sample. 

This figure also confirms our suspicion that it is the scale radius 
r_2 or concentration c v i r = r v i r /r_2 rather than the slope a that 
varies with n. We find that haloes forming in the n — —0.5 model 
tend to be more concentrated (cf . figure|4]below) than haloes forming 
in runs with steeper spectral indices. Therefore fits with a generalised 
NFW profile (cf. equation[9]l tend to favour smaller values of a for 
steeper n because these haloes tend to be less concentrated, and 
so we resolve the profile to smaller fractions of r_2, where the 
flattening of the profile is more apparent. Therefore a shallower 
effective slope a will tend to be preferred. 

It is a relatively straightforward exercise to obtain expressions 
for equation [TT] for the NFW profile and the |Moore et al.| < |1998) 
profile. Two other analytical model profiles have been promisingly 
applied to halo density profiles, the Eina sto] ( |1965| > and Prugniel 
|& Simien|([l 997) profiles, which provide better fits than the NFW 
profile. |Merritt et al7| q2Q06 ) argue that the Einasto model performed 
best in fitting halo profiles, followed closely by the Prugniel-Simien 
model. Expressions for 7 for the Einasto and Prugniel-Simien mod- 
els are given in the Appendix [B] 



We have taken great care to ensure that our estimates of r_2 are reliable. 
As haloes become more concentrated, r_2 can rapidly approach the inner- 
most reliably resolved radius, which can affect the accuracy with which r_2 
is estimated. This issue is particularly acute for haloes in the shallow-n cos- 
mologies. Therefore we do not consider haloes less massive than 0.75Af, . 
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Figure 3. The maximal inner slope as given by equation^^as a function of 
the radius normalised to r_2 for haloes in the M* sample. Upper panel: We 
plot all curves from all the simulations and over-plot a binned curve for each 
ran. Lower panel: Here only the averages of each simulation are shown and 
compared to a set of theoretical predictions derived from the NFW, Moore 
et al., Einasto and Prugniel-Simien profiles. Note that in this instance n 
corresponds to the shape parameter for the Einasto and Prugniel-Simien 
profiles. 



In the lower panel of figure [3] we over-plot the averaged 
7-curves with the theoretical predictions derived for the four 
analytic profiles mentioned above. We find that the Moore profile is 
unable to reproduce the observed behaviour. The NFW profile is 
consistent with our data for the 512-0.50 and 512-1.50 runs, but 
it fails to reproduce the continual flattening of 7 to small radii. In 
contrast, both the Einasto and Prugniel-Simien profiles capture the 
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models with spectral indices of n — —0.5, —1.0, —1.5, but who 
derived their concentrations from fits of NFW profiles. Although 
our concentrations are calculated in a non-parametric manner, it is 
reassuring that we see a similar trend to that reported in Navarro 
|et al.| (T997). The variation of concentration with mass (upper panel) 
suggests that this relation may be steeper in cosmologies with shal- 
lower spectral indices, which in turn may reflect the importance of 
merging for the growth of the most massive haloes in these cos- 
mologies. Although the trends are tentative, we believe that this 
is an extremely interesting figure and one that may allow to better 
understand the physical origin of the density profile. We shall return 
to this issue, albeit briefly, in §|4] 

The trend for concentration to decrease with steeper spectral 
index is more robust, and we show this explicitly by plotting the 
mean concentration of haloes in the M» sample versus spectral 
index (lower panel). This makes clear our assertion that it is the 
concentration that depends strongly on the spectral index. 
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Figure 4. In the upper panel we show the concentration of all haloes with 
mass larger than M* as a function of the number of particles within the 
virial radius. Additionally we plot a binned distribution of the concentration. 
The lower panel gives the mean concentration for objects in our Af* sample 
depending on the spectral index n. 

behaviour of our data well at small radij^j Interestingly we note that 
all of the analytical profiles tend to overestimate the slope of the 
density profile at large r/r_2, which appears to roll over and flatten 
off. This is most apparent for the data points from the n = —0.5 run. 

In Figure|4]we make explicit the connection between r_2 and 
the concentration c v ir = r v ir/V-2, showing how c v h- varies with 
halo mass (given by the number of particles within the virial radius 
iV v ir; upper panel), and the spectral index n (lower panel). A similar 
figure can be found in Navarro et al. ( 1997 ), who looked at scale-free 

5 Note that in this instance n corresponds to the shape parameter for the 
Einasto and Prugniel-Simien profiles. 



4 SUMMARY & CONCLUSIONS 

The primary motivation for this paper has been to establish what 
dependence, if any, the logarithmic slope a of the inner mass density 
profile of dark matter haloes has on the spectral index n of the linear 
matter power spectrum. For this purpose we have run a series of high 
resolution cosmological TV-body simulations of scale-free models 
(i.e. P(k) ~ Ak") with values of the spectral index n varying be- 
tween —0.5 ^ n ^ —2.75. By using scale-free models and fixing 
n in this manner, the problem of identifying correlations between 
a and n becomes a relatively straightforward one. Relatively, be- 
cause there is some freedom in the choice of criteria one can use to 
determine when to start and when to finish scale-free simulations. 

To address this, we have derived clear, reproducible and physi- 
cally motivated criterion that allows us to set up scale-free simula- 
tions and that is presented in § |2.1| In short, we start the simulations 
with the same integral power in the fc-range defined by the the total 
number of particles TV; it is especially important to have an objec- 
tive criterion for starting simulations when studying mass profiles. 
To determine when to stop our simulations, we follow [Navarro| 
et al. (19971 and use the evolution of the typical collapsing mass 
M*. We do so once the mass in a typical M* halo corresponds to 
approximately 42000 particles. 

Having established a set of well-defined criteria to set up and 
run cosmological simulations of scale-free models, we performed a 
sequence of high resolution runs (l/i _1 Mpc boxes, 512 3 particles) 
that we used to investigate the dependence of the central logarithmic 
slope of the dark matter halo density profile a on the spectral index 
n. We varied n between —0.5 and —2.75 and selected samples of 
well resolved haloes (N v i T ^ 3.15 x 10 4 ) in dynamical equilibrium 
in each run (ranging from ~ 270 haloes in the n — —0.5 run to 
~ 40 haloes in the n = —2.75 run). Using \ 2 fits to a generalised 
NFW profile, we identified preferred values of the inner slope a and 
indeed, we found a trend for the inner slope to become shallower 
with steeper spectral index - from a — 1.6 for n = —0.5 to 
a ~ 1.2 for n = -2.75. 

However, we argue that it is not the central slope a that 
depends on spectral index; rather, it is the scale radius of the 
halo, r s , or our preferred measure r_2, the radius at which the 
differential mass profile pr 2 reaches its maximum value. We 
have shown that haloes in different n models have similar radial 
profiles of the maximum slope 7 = 3(1 — p/p) when normalised 
to r_2. However, as already shown by Navarro e t al.| ( |1997[ l, 
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haloes that form in models with steeper n tend to be less centrally 
concentrated than haloes forming in models with shallower n, 
and so their mass profiles can be resolved to smaller fractions of 
r_2, where the flattening of the profile is more apparent. Haloes 
in these models would then appear to have shallower central profiles. 

We noted in the introduction that using scale-free simulations 
to study the effect of the spectral index of the power spectrum on 
the central structure of dark matter haloes was preferable to the 
approach t aken by |Ricotti| (|2003| >, |Cen et al.| ( |2004| >, |Colm et al.| 
(2004 ), and Ricotti et al. (2007 1, in which the box size and analysis 
redshift were varied to capture the behaviour of different spectral 
indices. The claims of Ricotti et al. (2007 ) are of particular interest; 
these authors argue that halo concentration is a universal constant 
and that dwarf galaxies identified at z=10 have logarithmic slopes 
shallower than 0.5. 

Our results strongly disagree with these claims. Although the 
central slope a that we obtain by fitting a generalised NFW profile 
does vary with n, this does not imply that the shape of the profile 
is sensitive to n, for the reasons presented above. Moreover, as 
we show in Figure[3] density profiles normalised by r_2, which 
is equivalent to concentration, have similar shapes on average, 
independent of spectral index; this implies that it is r_2 and 
consequently concentration that depends on n. It is possible that our 
profiles for y(r) may roll over and approach different asymptotic 
value at smaller radii than we can resolve with our simulations, 
but over the range of radii that we can reliably resolve we find no 
evidence for such behaviour. It should be noted that we checked for 
this by producing figure[3]for haloes in the high-mass sample (cf. 
Table [2} but we did not find any evidence for convergence. 



What can we learn about the physical origin of the density 
profile from our simulations? As we noted above, haloes in shallow - 
n models grow predominantly by accretion, whereas those in steep- 
n models grow via merging. This link between spectral index and 
accretion rate can be made explicit (see, for example, |Lu et al.] 
|2006) >, which makes the scale-free simulation a very interesting 
tool with which to study the origin of the profile. The results - 
that concentrations tend to be higher in shallow-n models, and 
that concentration decreases with increasing virial mass in these 
models (assuming that merging plays a role in the growth of the 
most massive halos) - indicate that "two-phase accretion" may be 
very important. However, to properly address this question a more 
detailed study is required. We are currently undertaking such a study, 
an account of which will be presented in a future paper. 
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We briefly consider the dependence of r_2 on n. We show in 
figure [4] that haloes forming in cosmologies with steeper spectral 
indices tend to be less centrally concentrated than haloes forming in 
cosmologies with shallower spectral indices, confirming and extend- 
ing results presented by Navarr o et al.|(T997) . This is interesting 
because we expect that haloes forming in shallow-n models do 
so predominantly by accretion, whereas haloes in steep-n models 
form via merging. If we examine the mass concentration relation 
in shallow-n models (upper panel in figure |4](, we find that more 
massive haloes tend to be less centrally concentrated than their less 
massive counterparts. In contrast, we find that this relation is far less 
pronounced in steep-n models. 

Prescriptions for halo concentration, such as Ek e et al.j ( [2001) , 
assert that concentration is determined by the halo's collapse redshift 
- haloes that collapse earlier do so when the mean density of the Uni- 
verse is higher, and so they tend to be more centrally concentrated. It 
has been argued that this reflects the rate at which the proto-halo ac- 
cretes mass (cf. Zhao et al. 2003 1. Haloes are said to assemble their 
mass in two distinct phases - fast-accretion phase during which the 
depth of the potential well is set, and a slow-accretion phase during 
which mass is added to the outer parts of the halo (e.g |Lu et al.|2006| l. 
This has important implications for the concentration. For example, 
|Wechsleret al. (2002) argued that during the fast accretion phase the 
central density is set by the background density, but as the accretion 
rate slows and the halo enters its slow accretion phase, the central 
density stays approximately constant and the halo's concentration 
grows in step with the virial radius. Merging also plays an important 
role; for example, Manriq ue et al.| ( |2003| > looked at the effect of 
major mergers on concentration and found that violent relaxation 
drives the remnant's density profile towards a form that is essentially 
identical to the one it would have formed through pure accretion. 



This paper has been typeset from a Tr-;X/ ETr-X file prepared by the 
author. 
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APPENDIX A: MAXIMAL EXPANSION 

A natural stopping criterion for a scale-free cosmological simulation 
is given by the requirement that the fundamental mode (or a multiple 
of it) turns non-linear; this defines a maximal expansion a max (cf. 
Efstat hTou et al.|1988) . But we like to remind the reader that in 
this study we compared simulations run with varying n values and 
hence had to define a different stopping criterion that allows us to 
cross-compare those runs at a similar evolutionary stage. 

However, in order to determine the maximal allowed expansion 
to obtain the fiducal M*^ value we follow Efstathiou et al.|p988] l 
and define a critical wavenumber ko, corresponding to a filtering 
scale for which the rms density fluctuations approach unity 



Table Al. The maximal allowed expansion factors. 
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This equation divides the power spectrum in two parts: for all fre- 
quencies below the critical frequency fcn the growth is still linear 
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and for all higher frequencies the growth is non-linear. Solving 
equation |Al| for ko yields: 



k = k 



Nv 



— Aa 1 
4 



n + 3 
2 



(A2) 



Therefore the evolution of the critical frequency ko with the expan- 
sion factor a of the universe is now known and can be inverted to 
calculate a maximal a max as a function of a pre-defined threshold 
wavenumber feo 



-Ax" 



n + 3 



where x is given by the relation 

ko = xk-^y 



(A3) 



(A4) 



As the threshold frequency we choose the smallest wavenumber 
that we are able to resolve numerically, i.e. the fundamental mode 
in our computational domain as given by x = 2/N . The maximal 
allowed expansion factors for our scale-free models are listed in 
Table | A 1 1 and we encourage the reader to compare them against 
the adopted expansion factors summarized in Table [T] Our stopping 
criterion complies with the requirement that the fundamental mode 
is still in the regime of linear growth. 



APPENDIX B: MAXIMUM SLOPE : EXPRESSIONS FOR 
EINASTO AND PRUGNIEL-SIMIEN MODELS 

For the Einasto model, the density profile can be written as 

PEin(r) = exp{-2n(f« (Bl) 

where p and f are the density and the radius in units of p_2 and r_2, 
respectively, and n gives the shape of the density profile. We will 
use this notation throughout this appendix. For formulations of the 
profiles in various units and conversions between them see |Merritt| 
|etllT| ( [2006) and |Graham eTaiJj20Q6l 

From this, the enclosed mass MEin can be obtained. This is 
given by: 

M Ein = 27re 2n (2n) 1_3n 7 (3n, 2nf » ) (B2) 
where 7(0, x) is the lower incomplete Gamma Function given by 

f{a,x) = f t a_1 e~'d£. (B3) 
Jo 



This yields 



7(f) 



= 3 



3 7 (3n,2nfi) 



(B4) 



for the maximum slope. 

For the Prugniel-Simien model, the density profile is written as 



Dark Matter Halo Profiles in Scale-Free Cosmologies 

p(r) = r~ p exp {-n{2 - p) (r™ - l)} (B5) 
where p is a function of n which can be approximated by 
p(n) = 1 - 0.6097/n + 0.05463/n 2 . (B6) 
The enclosed mass profile can be written as 

M (f ) = 7^ L p e n(2 ~ p) (n(2 - p)) 1 -^ 3 ^ 7 (n(3 - p), (2 - p)nf - ) (B7) 
which leads to 

7ffl=3 (,. ^2-^-^^ (B8) 
y 6 fP- 3 7 (n(3 -p), (2 -p)nfs J y 

for the maximum slope. 



